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Abstract. The paper presents a methodology for the enhanced stiffness analysis of parallel manip- 
ulators with internal preloading in passive joints. It also takes into account influence of the external 
loading and allows computing both the non-linear "load-deflection" relation and the stiffness ma- 
trices for any given location of the end-platform or actuating drives. Using this methodology, it is 
proposed the kinetostatic control algorithm that allows to improve accuracy of the classical kine- 
matic control and to compensate position errors caused by elastic deformations in links/joints due 
to the external/internal loading. The results are illustrated by an example that deals with a par- 
allel manipulator of the Orthoglide family where the internal preloading allows to eliminate the 
undesired buckling phenomena and to improve the stiffness in the neighborhood of its kinematic 
singularities. 
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1 Introduction 

Parallel manipulators have become very popular in many industrial applications 
due to their inherent advantages of providing better accuracy, lower mass/inertia 
properties, and higher structural rigidity compared to their serial counterparts (TJ. 
These features are induced by the specific kinematic structure, which eliminates 
the cantilever-type loading and allows to minimize deflections caused by external 
torques and forces. One recent development in this area, which is targeted at high- 
precision manipulation, is a replacing the standard passive joints by preloaded ones, 
which contain internal passive springs eliminating the backlash or ensure some de- 
gree of static balancing [2 3 J. This modification obviously improves the manipula- 
tor performances but requires some revision of existing stiffness analysis techniques 
that are in the focus of this paper. 

In most of previous works, the manipulator stiffness analysis was based on the 
linear modeling assumptions which ignore influence of the external or internal 
forces 10] |5] |6] |7] ID . Consequently, relevant techniques are targeted at linearization 
of the "force-deflection" relation in the neighborhood of the non-loaded equilibrium, 
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which is perfectly described by the stiffness matrix J9][T0]. However, in the case of 
non-negligible internal and/or external loading, the manipulator may demonstrate 
essentially non-linear behaviour, which is not exposed in the unloaded case ifTTI . 
In particular, the loading may potentially lead to multiple equilibriums, to bifurca- 
tions of the equilibriums or to static instability of certain manipulator configurations 

MM- 

This paper presents an extension of our previous results [14] devoted to the 
stiffness analysis of parallel manipulators by generalizing them for case of inter- 
nal preloading [15| in the passive joints. It implements the virtual joint method 
(VJM) of Salisbary [16| and Gosselin [17| that describes the compliance of the 
manipulator elements by a set of localized multi-dimensional springs separated by 
rigid links and perfect joints. The proposed technique allows computing the loaded 
equilibrium, finding the full-scale "load-deflection" relation and evaluating the cor- 
responding stiffness matrices for any given location of the end-platform or actuating 
drives ifTHl . It is also developed a kinetostatic control algorithm that allows to im- 
prove accuracy of the classical kinematic control and to compensate position errors 
caused by elastic deformations in links/joints due to the external/internal loading. 

The remainder of this paper is organized as follows. Section 2 defines the re- 
search problem and basic assumptions. Section 3 deals with computing of the loaded 
static equilibrium and corresponding "load-deflection" relation. Section 4 focuses 
on its linearization and evaluation of the stiffness matrix. Section 5 presents the 
kinetostatic control algorithm. Section 6 contains an illustrative example. And fi- 
nally, Section 7 summarizes the main results and contributions. 



2 Manipulator model 

Let us consider a general parallel manipulator that is composed of n serial kine- 
matic chains connecting a fixed base and a moving platform Figure 1 . It is assumed, 
that the chain architecture ensures kinematic control of the manipulator but may 
introduce some redundant constraints that improve the rigidity. Following the VJM- 
concept [17 1, let us presents the manipulator chains as sequences of pseudo-rigid 
links separated by rotational or translational joints of one of the following types: 
(i) perfect passive joints ; (ii) preloaded passive joints that include auxiliary flexible 
elements; (iii) virtual flexible joints that describe compliance of the actuators and 
manipulator links; (iv) actuating joints. Using this notation the geometrical model 
of the chain may be written as 

t = g (p,q,tf,e), (1) 

where the vector t = (p, (f>) r includes the Cartesian position p = (x, y : z) and 
orientation <p = (<p v , <p y , (p z ) of the end-platform, p is the vector of actuated coor- 
dinates (they are constant for static analysis), the vector q contains coordinates of 
all perfect passive joints, the vector # includes coordinates of the preloaded pas- 
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Fig. 1 Typical parallel manipulator and VJM-model of its kinematic chain. 
(Ac - actuator; Ps - Passive joint) 




(a) Perfect passive joint (b) Passive joint (c) Passive joint 

with a linear spring with a non-linear spring 



Fig. 2 Examples of auxiliary springs in preloaded passive joints. 



sive joints, and the vector 9 collects coordinates of all virtual springs describing 
elasticity of the links and joints. 

The above mentioned elements of the kinematic chain differ in their static char- 
acteristics. In particular, the joints (i) and (iii) are described by the standard expres- 
sions lfl4l 



where x q and Xg are the generalized force/torque reactions corresponding to the 
aggregated vectors of the passive joint coordinates q and virtual joint coordinates ; 
Kg is the generalized stiffness matrix of all virtual springs. However, the preloaded 
passive joints (ii) may include both linear and non-linear auxiliary springs, some 
examples of which are shown in Figure [2] In this paper, we will describe statics of 
the preloaded joints by a general expression 



where T$ is the generalized force/torque reactions corresponding to the aggregated 
vectors of the preloaded joint coordinates i9n defines the preloading value; K# 
is the generalized stiffness matrix of preloaded joints, and the vector function h(...) 
is assumed to be piecewise-linear, such that each of its scalar components 
can be expresses either as the difference (#,■ — #,o), or its positive or negative part 
— #;o] + , [#/ — #;0p ( see Figure 2 for details). 

Using these assumptions, let us derive the stiffness model of the considered ma- 
nipulator and sequentially consider the following sub-problems: (i) computing the 
loaded static equilibrium and obtaining the "load-deflection" relation; (ii) lineariza- 



i q = and Tg = Kg • 9 



(2) 



t#=K# •h(#-i9 ) 



(3) 
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tion of this relations in the neighborhood of this equilibrium and computing the 
stiffness matrix; (iii) developing the kinetostatic control algorithm, which allows to 
compensate position errors caused by the elastic deformations and preloading. 



3 Static equilibrium 

Let us obtain first the configuration of each kinematic chain (q, 9, ■&) and external 
force F that correspond to the static equilibrium with the end-point location t. Ob- 
viously, it is a dual problem compared to the classical static analysis but it is more 
reasonable here because of strictly parallel structure of the considered manipulator 
(see Figure 1). The latter allows applying the same technique to all kinematic chains 
(with the same end-point location) and to compute the total external loading as the 
sum of the partial loadings. 

Taking into account the assumption on the piecewise-linear property of the func- 
tion h(.), let us perform regrouping of the variables. In particular, for each current 
configuration of the chain, the coordinates of the preloaded passive joints described 
by the vector # may be separated into two parts # e and # ? , where the first one 
corresponds to the active state of the auxiliary springs and the second part describes 
non-active springs (see Figure 2 for geometrical interpretation). This allows replac- 
ing the original set of the configuration variables (q, 6, ■&) by a set of two vectors 
(q, 9), where q aggregates the joint coordinates (q, # 9 ) that currently are passive 
and the vector 9 collects all spring coordinates (9, -&q) (both virtual and passive). 

Using these notations and applying the virtual work technique, the static equilib- 
rium equation of the kinematic chain may be written as 

r e -F = K e -(e-e y, j£-f = o (4) 

where F is the external force applied at the end-point of the chain, the vector 9q = 
[0 T , Pq] aggregates the spring preloadings (which is obviously zero for the virtual 
springs), Kg = diag(Kg,K^), and Jq, J q are the kinematic Jacobians derived from 
(1) by differentiating it with respect to 9, q. This system of equation (4) combined 
with the geometrical model (1), which must be rewritten in terms of the redefined 
variables 

t = g(q,0). (5) 

This yields the desired joint coordinates of the static equilibrium for a separate kine- 
matic chain with given end-point location. 

Since the derived system is highly nonlinear, in general case a desired solution 
can be obtained only numerically. In this paper, it is proposed to use the following 
iterative scheme 
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where the starting point (6>o,qo) is also computed iteratively, started from a near- 
est unloaded configuration where the joint coordinates are easily obtained from the 
inverse kinematic model. On the following iterations, to improve convergence, the 
system variables are slightly randomly disturbed. As follows from computational 
experiments, the proposed iterative algorithm possesses rather good convergence 
(3-5 iterations are usually enough). 



4 Stiffness matrix 

To compute the desired stiffness matrix, let us consider the neighborhood of the 
equilibrium configuration and assume that the external force and the end-effector 
location are incremented by some small values 8F, 8t. Besides, let us assume that 
a new configuration also satisfies the equilibrium conditions. Hence, it is neces- 
sary to consider simultaneously two equilibriums corresponding to the manipulator 
state variables (F,q,0,t) and (F + <5F,q + Sq, 9 + 59, t + <5t). Relevant equations 
of statics may be written as 



j£F = K e (0-0 o ); J^F = 0; 

(Je + 8Je) T (F + 5F) = K e (9-9 Q + 86) ; (j q + 5J„) r (F + <5F) = 

where 8j g (q, 9) and 5Je(q, 9) are the differentials of the Jacobians due to changes 
in (q, 9). Besides, in the neighborhood of (q, 9), the kinematic equation (5) may be 
also presented in the linearized form: 

5t = J e (q, 9)-89+J q (q, 0)-8q, (8) 
Hence, after neglecting the high-order small terms and expanding the differen- 

~ T 

tials via the Hessians of the function f = g(q, 9) F 



H F qq = d 2 ¥/dq 2 ; n F ee = d 2 V/d~9\ U F e = (U F 9q ) T = d 2 ¥/dq dd, (9) 
equations (7) may be rewritten as 

J T e (q, 9)-8F + H F q (q, 9) ■ 8q + H F e (q, 9)-89=K e - 89 



(10) 

J^(q, 9) ■ SF + n F qg (q, 9)-8q + H F e (q, 6) ■ 86 = 
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Besides, here the variable 89 can be eliminated analytically: 89 



H-5F- 



k£ • Hq 9 • Sq, where k£ = (K e 
with unknowns 5F and 8q 



Ho 



. This leads to a system of matrix equations 
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(11) 



from which the desired Cartesian stiffness matrix of the chain K c may be obtained 
by direct inversion of the the left-hand side and extracting from it the upper-left 
sub-matrix of size 6x6: 







* * 





(12) 



Je -k ■ J Jq + Je 
J q +H 9 g -kg • J e ¥L qq + fl q9 -k e 

Finally, when the stiffness matrices for all kinematic chains are computed, the 
stiffness of the entire multi-chain manipulator can be found by simple summation 
Kj; = Y."i=\ Krf- It should be noted that, because of presence of the passive joints, the 
stiffness matrix of a separate serial kinematic chain is always singular, but aggrega- 
tion of all the manipulator chains of a parallel manipulator produce a non-singular 
stiffness matrix. 



5 Kinetostatic control 

In robotics, the manipulator motions are usually generated using the inverse kine- 
matic model that allows computing the input (reference) signals for actuators p cor- 
responding to the desired end-effector location t. However, for manipulators with 
preloaded passive joints, the kinematic control becomes non-applicable because of 
changes in the end-platform location due to the internal loading. Hence, in this case, 
the control must be based on the inverse kinetostatic model that takes into account 
both the manipulator geometry and elastic properties of its links and joints lfl2l . 

Using results from the previous sections, the desired inverse kinetostatic trans- 
formation can be performed iteratively, in the following way: 

Step#l. For given target location of the end-platform t, compute initial values of 
the actuated coordinates po by applying the inverse kinematic transformation. 

Step#2. For current values of the actuated coordinates p ( - and target location of 
the end-platform t, find the equilibrium configuration for each kinematic chain 
and compute the corresponding total external loading required to achieve the 
target location. 

Step#3. If the computed external loading is less than the prescribed error, i.e. 

Fj < Ef , stop the algorithm, otherwise continue the next step 
Step#4. Repeat Step#2 several times in the neighborhood of the current solution 

Pi and evaluate numerically the matrix S' Fp = dF' E /dpj describing the sensitivity 

of F with respect to p . 
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(a) 3D Orthoglide (a) 2D Orthoglide 

Fig. 3 Architecture of the Orthoglide manipulator and its planar version. 

Step#5. Compute new value of the actuated coordinates p, + i = p,- — S' Fp 1 • 
and repeat the algorithm starting from Step#2. 

As follows from simulation results, this algorithm demonstrates good conver- 
gence and can be used both for on-line and off-line trajectory planning. It was suc- 
cessfully applied to the case-study presented in the following Section. 



6 Application example 

Let us apply the proposed techniques to the stiffness analysis of the planar manip- 
ulator of the Orthoglide family (Figure 3). For illustration purposes, let us assume 
that the only source of the manipulator elasticity is concentrated in actuated drives, 
while the passive joints may be preloaded by (i) standard linear springs, or (ii) non- 
linear springs with mechanical stop-limit (see Figure 2 for details). 

For this manipulator, the kinematic model includes a single parameter L (the 
leg length) and the dexterous workspace was defined as the maximum square area 
that provides the velocity (and force) transmission factors in the range [0.5, 2.0]. 
Using the critical point technique developed for this type of manipulators 1 1 9 ] , it was 
proved that the desired square vertices are located in the points Qi(—p, — p) and 
Qi(p, p), where p = 0.45 L. Besides, the square centre Qo(0, 0) is isotropic with 
respect to the velocity and force transmission. The parameters of the actuating drives 
are also assumed identical and their linear stiffness is denoted as Kg . The auxiliary 
springs incorporated in the passive joints adjacent to the actuators are described 
by two parameters: the angular stiffness coefficient K$ and the activation angle #o 
that defines the preloading activation point. During simulation, the manipulator end- 
point was displaced by value A in the direction QoQi or Q0Q2, and it was computed 
corresponding magnitude for external force F. 

The stiffness analysis results are summarized in Figures 4, 5 and in Table 1 . As 
follows from them, the original manipulator (without preloading in passive joints) 
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Fig. 4 Force-deflection relations F = f(A/L) in critical points: 
(l)^ tf =0; (2) K$ = 0.01 Kg L 2 ; (3) K# = 0. 1 K e L 2 
(case of preloading with linear springs). 
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Fig. 5 Compliance maps for cases of: (a) manipulator without preloading; 
(b) manipulator with preloading non-linear springs with K$ = 0.5 Kg L 2 and t9rj = Jt/12 (b). 



demonstrates rather low stiffness in the neighborhood of the point Q2, which is 
roughly 4 times lower than in the isotropic point Qq. In contrast, the linear stiffness 
in the point Q\ is twice higher than in the point Qq. Besides, in the point Q2, the 
external loading may provoke the buckling phenomenon that is caused by a local 
minimum of the force-deflection relation. In this case, the distance-to-singularity is 
essentially lower that it is estimated from the kinematical model and the manipulator 
may easily loose its structural stability. 

To improve the manipulator stiffness and to avoid the buckling in the neigh- 
borhood of Q2, the passive joints were first preloaded by linear springs with ac- 
tivation angle #0 = 0. As follows from Figure 4, the preloading with parameter 
K$ = 0.1 Kg L 2 allows completely eliminate buckling and improves the stiffness by 
the factor of 2.3. On the other hand, the stiffness in the points Qq and Q\ changes 
non-essentially, by 10% and 5% respectively. Hence, with respect to the stiffness, 
such preloading has positive impact. 

The only negative consequence of such preloading is related to changes of the 
actuator control strategy. In fact, instead of standard kinematic control, it is nec- 
essary to apply the kinetostatic control algorithm presented in Section 5. It allows 
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Table 1 Manipulator stiffness for different linear preloading. 



Stiffness in preloaded joints K$ = 0.01 Kg L? 


0.05 Kg L 2 


0.1 Kg L 2 


Point go (isotropic point) 


Actuating joint coordinates p L L 
Manipulator stiffness K £ . Kg 1.01 Kg 


L 

1.05 Kg 


L 

1.10 Kg 


Point Q\ (neighborhood of "bar" singularity) 


Actuating joint coordinates p 0.437 L 0.433 L 
Manipulator stiffness K c 2.276 Kg 2.286 Kg 


0.419 L 
2.329 Kg 


0.402 L 
2.382 Kg 


Point Qi (neighborhood of "flat" sit 


igularity) 




Actuating joint coordinates p 1.345 L 1.356 L 
Manipulator stiffness K c 0.24 Kg 0.27 Kg 
Critical force F cr 0.020 Kg L 0.027 Kg L 


1.399 L 
0.39 Kg 


1.453 L 
0.55 Kg 



compensating the position errors caused by elastic deformations due to the internal 
preloading and to achieve the target end-point location with modified values of the 
actuated joint coordinates. As follows from Table 1, corresponding adjustments of 
the joint coordinates may reach 0. 1 L and are not negligible for most of applications. 

The most efficient solution that eliminates this problem is using of non-linear 
springs with mechanical stop-limits that are activated while approaching to Qi. For 
instance, as follows from dedicated study, the preloading with the parameters K# = 
0.5 Kg L 2 , i?o = ft/12 provides almost the same improvements in Q2 as the linear 
spring while preserving usual control strategies if the preloading is not activated. 
The efficiency of this approach is illustrated by the compliance maps presented in 
Figure 5. 



7 Conclusions 

Recent advances in mechanical design of robotic manipulators lead to new parallel 
architectures that incorporates internal preloading in passive joints allowing to im- 
prove accuracy but leading to revision of existing stiffness analysis techniques. This 
paper presents new results in this area that allow simultaneously evaluate influence 
of internal and external loading and compute both the non-linear "load-deflection" 
relation and the stiffness matrices for any given location of the end-platform or 
actuating drives. Using this methodology, it is proposed the kinetostatic control al- 
gorithm that allows to improve accuracy of the classical kinematic control and to 
compensate position errors caused by elastic deformations in links/joints due to the 
external/internal loading. The efficiency of this technique is confirmed by an appli- 
cation example that deals with a parallel manipulator of the Orthoglide family where 
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the internal preloading allows to eliminate the undesired buckling phenomena and 
to improve the stiffness in the neighborhood of its kinematic singularities. 

In future, these results will be generalized to other types of preloading that may 
be generated by external gravity-compensation mechanisms and also applied to mi- 
cromanipulators with flexure joints. 
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